Nonlinear modes in the harmonic T'T-symmetric potential 
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We study the families of nonlinear modes described by the nonlinear Schrodinger equation with 
the 'PT-symmetric harmonic potential x 2 — 2iax. The found nonlinear modes display a number of 
interesting features. In particular, we have observed that the modes, bifurcating from the different 
eigenstates of the underlying linear problem, can actually belong to the same family of nonlinear 
modes. We also show that by proper adjustment of the coefficient a it is possible to enhance 
stability of small-amplitude and strongly nonlinear modes comparing to the well-studied case of the 
real harmonic potential. 
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I. INTRODUCTION 

The interest in the stationary modes of the nonlinear 
Schrodinger equation with a potential has been raised 
about two decades ago in connection with the applica- 
tions to the meanfield dynamics of Bose-Einstein con- 
densates pj and later on in the context of in the context 
of optical applications @ and in particular of propaga- 
tion of dispersion managed solitons in fibers Q. Various 
aspect of the nonlinear modes in a parabolic trapping 
potential with homogeneous 0-0] and inhomogeneous [|[ 
nonlinearities have been intensively studied. A compre- 
hensive analysis of the structure of the nonlinear modes 
and their stability can be found in @43- Further, ac- 
counting that interaction of a particle with a potential 
in practice is not absolutely elastic, and energy losses 
are possible, in [i| there have been addressed nonlinear 
modes in a complex parabolic potential (1 — i)x 2 sup- 
ported by a homogeneous gain. Due to its dissipative 
nature, the complex parabolic potential has properties 
very different comparing to its real counterpart. In par- 
ticular, for the fixed parameters of the dissipative model, 
the stable nonlinear modes appear as isolated attractors 
and do not constitute continuous families. Another inter- 
esting feature of the complex parabolic potential is that 
in the limit of the strong defocusing (or repulsive) non- 
linearity, the so-called Thomas-Fermi approximation of 
the model is described by the balance between the losses 
and the gain. This is not the case of the conservative 
potential, where the behavior of the nonlinear modes in 
the Thomas-Fermi limit is determined by the balance be- 
tween the dispersion (or diffraction) and the nonlinearity. 

In the meantime, recently there has appeared a rapidly 
increasing interest [loj in linear and nonlinear proper- 
ties of the systems with potentials obeying the so-called 
PT symmetry. This interest was initiated by the pa- 
per , and more recently by the experimental observa- 
tion of PT symmetry breaking in optics [12J, as well as 
by several theoretical suggestions of realization of "PT- 
symmetric optical systems 13 1. 



The nonlinear extensions of the PT-symmetric struc- 
tures have been first considered in Later on the 
nonlinear modes have been studied in the periodic [l5j . 
Gaussian [l6[, and sech 2 -shaped [TtJ PT-symmetric po- 
tentials, as well as in the harmonic trap with rapidly 
decaying PT-symmetric imaginary component (l8| . We 
also mention studies of gap solitons in PT-symmetric 
optical lattices combined with real superlattices [l9[ and 
optical defect modes in PT-symmetric potentials [20|. 
The modes and their stability in the systems with PT- 
symmetrically modulated nonlinearity landscapes have 
been recently reported in 



21] 



It turns out, however, that nonlinear modes in the PT- 
symmetric parabolic trap have not received any atten- 
tion, so far, while such a potential, namely [x — ia) 2 , 
has been introduced and well studied in the linear the- 
ory [22|-[24j]. Meantime, as it will be shown below, the 
nonlinear modes in the PT-symmetric harmonic poten- 
tial display rather unusual properties, which cannot be 
observed either in conservative or in dissipative poten- 
tials of a general kind. The main goal aim of the present 
work is to perform a detail study of such modes. 

The rest of the paper is organized as follows. In the 
next section we introduce the nonlinear model with the 
PT-symmetric harmonic potential and briefly discuss its 
physical relevance. In Sec. Mil we discuss some properties 
of the underlying linear model. Next, in Sees. HVl andlVl 
and we report the families of nonlinear modes, as well as 
a detail investigation of their stability. Sec. IVII concludes 
the paper. 



II. THE MAIN MODEL 

Our main object in the paper is the nonlinear 
Schrodinger equation with a PT-symmetric parabolic po- 
tential: 



KJz 



-q xx + (x 2 - 2iax)q - <r\q\ 2 q, 



(1) 
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where a > and a = 1 and a = — 1 correspond to focus- 
ing and defocusing nonlinearity (hereafter we use termi- 
nology relevant to optical applications). Physically the 
dimensionless Eq. (fT]) naturally appears as an equation 
modeling a beam guidance in a medium whose refractive 
index n(x) — n r (x) + irii(x) has parabolic modulation of 
the real part n r {x) = x 2 and linear modulations of the 
imaginary part rii{x) = lax. Even more generally, any 
smooth enough symmetric profile of the refractive index 
n r (x) = n r (—x) and anti-symmetric modulation of its 
imaginary part, n^x) — — rij(— x) leads to the model ([T]) 
if a guided beam is narrow enough allowing for the use 
of the first order terms of the Taylor expansion of the 
complex index n{x). 

In this paper we are interested in stationary modes, 
which are searched in the form q(z, x) — w(x)e l ^ z , where 
P is the propagation constant. We consider localized so- 
lution which obey the zero boundary conditions: 

lim \q{z,x)\ = 0. (2) 



For the next consideration it is convenient to introduce 
the representation j3 = b— a 2 , where & is a new parameter, 
which allows one to arrive at the following stationary 
equation: 

w X x — bw — (x — ia) 2 w + a\w\ 2 w = 0. (3) 

Recalling that the existence of the modes implies the 
balance between the diffraction and the nonlinearity, as 
well as between gain and losses, we also rewrite Eq. ^ 
in the hydrodynamical from 

Pxx - (b- a 2 - x 2 )p + <rp 3 - = 0, (4a) 

P A 

j x = -2axp 2 7 (4b) 

where p(x) = \w(x)\ is the field modulus, while j(x) — 
0x{x)p 2 (x), with 9(x) = argw(x), is the real-valued cur- 
rent. From (j4]) one readily concludes that both p(x) and 
j(x) — 6 x {x)p 2 {x) are even functions. The current j(x) 
has a local maxima at x = 0, while p(x) has either a 
local maximum or a local minimum at x = 0. Moreover, 
it follows from Eq. ([2]) that j — > at x — > oo, and hence 
taking into account that j x (x) does not change sign for 
x 0, we deduce from Eq. (|4a[) that j(x) does not be- 
come zero at any finite x, and hence the same is valid 
to p(x) [since otherwise the last term in Eq. (|4a|) would 
give a singularity] . The absence of zeros of the field con- 
trasts to the known behavior of the nonlinear modes in 
a real harmonic potential, while is known for the linear 
■pT-symmetric modes [231 ]. which are briefly outlined in 
the next section. 



III. LINEAR MODES 

Let us recall some relevant properties of the linear 
problem 

d 2 

L n W n = 0, C n = -J^, ~ bn ~ ( X ~ > ( 5 ) 

which can be formally obtained by setting a = in 
Eq. ([3]). Hereafter a tilde distinguishes solutions of the 
linear problem. The set of the eigenvalues of the problem 
([5]) does not depend on a and consists of an equidis- 
tant sequence b n = — (2n + 1), n = 0,1, — Cor- 
responding eigenfunctions can be written as w n (x) — 
Cn^nix — ia), where ^„(a;) = H n {x)e~ x I 2 is the nth 
Gauss-Hermite mode, J , 4' n {x)ijj% l {x)dx — S n , m y/ir2 n n\, 
H n (x) is the nth Hcrmite polynomial, and c„ are the 
positive coefficients providing the normalization condi- 
tion J w n (x)Wn(x)dx — 1 (hereafter we omit the inte- 
gration limits wherever the integration is over whole real 
axis, and the asterisk denotes complex-conjugation). 

Unlike in the conservative case a = 0, for a > 
the eigenfunctions w„(x) are not orthogonal. Us- 
ing the relation (see e.g. [Hj]): H n (x + x ) = 
El^qC^oT^Hkix) where = n\/[k\(n — k)\] are 
the binomial coefficients), for any n and m one finds 

/ w n (x)wl l (x)dx = c n c m e a y/irx 

ELo C*C; n 2 k k\(-l) n - k (2ia) n+m - 2k = 

c n c m e a2 v ^2 :l± ^ ±R i 3n + m a9 p lL { p 9 \~2a 2 ), 

wherep = min(n,m), g = \n— m\, and Lp (x) is the gen- 
eralized Laguerre polynomial. Setting n — m we obtain 
the expression for the normalization coefficients c n : 

-q 2 /2 

Cn = , (6) 

V0F2»n!L n (-2«2) 

For a = the eigenfunctions w n (x) are real- valued 
(up to irrelevant phase shift). Moreover, w n (x) is an 
even (odd) function if n is even (odd). For a / the 
eigenfunctions are complex-valued and are neither even 
nor odd. Instead, they can be chosen to have even real 
part and odd imaginary part. 

IV. BIFURCATIONS OF NONLINEAR MODES 

Turning now to the nonlinear problem, we observe that 
the eigenvalues b n , n = 0, 1, . . . are the bifurcation points 
where families of nonlinear modes branch off from the 
zero solution w(x) = 0. The nonlinear modes w n (x) be- 
longing to the nth family have the same symmetry as the 
corresponding linear eigenfunction w n (x) . In the vicinity 
of the n-th bifurcation point, the nonlinear modes w n (x) 
can be described by means of asymptotic expansions 

w n (x) = sw n + C(e 3 ), b n = b n + (re 2 b (2) + o(e 2 ), (7) 
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where £ < 1 is a formal small parameter. Since w n (x) 
were chosen normalized, in the leading order the total 
energy flow U — J \w n (x)\ 2 dx (hereafter all the integrals 
are taken over the whole real axis), is equal to e: U ~ e 2 . 
The solvability condition for the £ 3 -order equation yields 



, (2) = ! wl{x)w* n {x)dx 
Jw 2 n (x)dx 



(8) 



Since for a = the eigenfunctions w n {x) are real val- 

(2) 

ued, one has that 6„ ' > for any n. For a > the 
eigenfunctions w n {x) are complex- valued. However, par- 

(2) 

ity of their real and imaginary parts ensures that 6„ is 
nevertheless real for any n and a. It is straightforward 

(2) 

to obtain explicit expressions for b n . For the two lowest 
families (n = and n = 1) one has 



°0 — 



2tt 



6f > = 3 ' 



l + 2a 2 



1 + 2a 2 



(9) 



It follows from Eqs. (|9]) that b is positive for all a 

(2) 

while b\ is positive for small a, but becomes negative 
for a > \fl + \/2. Regarding the next families, we have 

(2) 

found that for n = 2 the coefficient b 2 changes sign 

(2) 

twice. For n = 3, however, the coefficient b y 3 ' changes 
sign only once, becoming negative for all sufficiently large 
a. 

From Eqs. © we also arrive at another interest- 
ing observation: 
tially fast with 
lim 



the coefficients b { grow exponcn- 



du I 



\b=bo,i 
1,(2) 



This, in particular, means that 
0. Taking into account that the 

coefficients b^' and &j have opposite signs for a 3> 1, 
one can expect that for large a the nonlinear modes bi- 
furcating from bo and b\ merge (or intersect) at some 
value of the energy flow U. 

The latter situation seems to be counterintuitive and 
strongly contrasting to what is known for the conser- 
vative harmonic potential, where the modes bifurcat- 
ing from different eigenstates of the linear problem do 
not merge. In order to check this issue we performed 
the direct numerical study of the families of nonlinear 
modes. The characteristic results are summarized in 
Fig. [TJ where the families of nonlinear modes are shown 
on the plane (6, all) for several different values of a. Re- 
spectively, the modes corresponding to the focusing (de- 
focusing) nonlinearity are situated above (below) the axis 
all — 0, which is indicated with the dashed line. 

For the sake of comparison, in the left upper panel 
of Fig. [T] we show the families of nonlinear modes for the 
well-studied real harmonic oscillator (3,[||,0]; which in our 
case corresponds to a = 0. Increasing a (see the other 
panels of Fig. Q]), we observe that already at a = 1 in the 
defocusing medium the nonlinear modes bifurcating from 
&o = —1 and b\ = —3 (as well as the ones bifurcating from 
&2 = — 5 and 63 = —7) indeed appear to be connected in 
a single family.For larger a (e.g. for a = 2) the structure 
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FIG. 1: The lowest families of nonlinear modes for different 
a. The fragments of curves corresponding to stable nonlinear 
modes are shown in bold. The nonlinear modes indicated with 
the points (a)-(d) in the panel a = 1 are explicitly shown in 

Fig.rj 
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FIG. 2: (Color online) The modulus p(x) and the current j(x) 
for stable nonlinear modes corresponding to a — 1. Panels 
(a)-(d) correspond to nonlinear modes indicated by the points 
(a)-(d) in the panel a = 1 of Fig. [T] 
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of the nonlinear modes becomes more complicated and 
the higher families (the ones bifurcating from 64 = —9 
and 65 = —11) also turn to be involved in creation of a 
singe family snaking through the linear eigenstates with 
n = 2, 3, 4, and 5. For a — 2 one can see the connection of 
the modes not only in the defocusing medium but also in 

the focusing one. Since a — 2 > \J 1 + y/2, Eqs. (|9]) imply 

that the coefficient b\ is negative for a = 2, and thus, in 
contrast to the cases a = 0,a = 0.15 and a = 1, the slope 
d(alf)/db is negative in the vicinity of the bifurcation 
from the point b\ . In Fig. [2] we show the field modulus 
p{x) and the superfluid current j{x) for several stable 
nonlinear modes corresponding to a = 1. In accordance 
with the discussion in Sec. [Ill both p(x) and j(x) are even 
functions, and for all the shown modes the current j{x) 
has a maximum at x = 0. The field modulus p(x) has a 
maximum at x — for the nonlinear modes (a), (c) and 
(d). For the nonlinear mode (b) the field modulus has a 
local minimum at x = 0. 

It is interesting to observe, that the described behav- 
ior of the modes allows one to suggest that it is possible 
to use continuous deformation to transform one of the 
modes of the conventional linear harmonic oscillator to 
another one having different parity. Indeed, to this end 
it is enough to properly change the strength of the non- 
conservative potential a and the intensity of the beam U. 
Notice that the stability of the modes, important for any 
practical realization of such a deformation is discussed in 
the next section. 



V. STABILITY OF THE NONLINEAR MODES 



Let us now analyze the spectrum of the operator L in 
the vicinity of the nth bifurcation point. In the linear 
limit (i.e. for e = 0) the operator L acquires the form 



£n 

-ct 



(11) 



where C n is defined in ([S]). The spectrum of the operator 
L„ consists of two sequences. Eigenvalues and eigenvec- 
tors of the first sequence read uj^\ = 2(n — k), p^\ = 
(wfe(cc),0) T , k = 0,1, — The second sequence reads 
-2(n-fc), p£? = (0,wt(x)) T , k = 0,1,.... 



u n,k 



A. Analytical results 



First, we notice that the operator L„ has a double zero 
eigenvalue u)n,n = <^n,n = 0. Generically, passing from 
the linear limit e = 0toe>0,a double eigenvalue 
splits into two simple eigenvalues. However, in the case 
at hand, the splitting of the double zero eigenvalue can 
not occur. Indeed, if the zero eigenvalue splits into two 
simple ones, they will be either both real and of opposite 
signs or complex conjugated. Either of these possibilities 
means that for e =/= the eigenvalue u) — is no longer in 
the spectrum of the operator L. This, however, contra- 
dicts to the established above properties of the operator 
L. Thus, for e 7^ the operator L also has the double 
zero eigenvalue. 

Besides of the double zero eigenvalue, the operator L„ 

has 2ra double eigenvalues: £l n ,k — u „ k ~ n 2n-fc' wnere 
k runs from to 2n except for k = n. Again, the double 
eigenvalue fl Ut k generically splits into two simple eigen- 
values, which will be either both real or complex con- 
jugated. At the same time, the opposite double eigen- 
value ri n ,2n-fc = — &n,k will split in the same manner. 
Since the double eigenvalues f2 ni /- and f2 nj 2n-fc behave 
in the same way, it is sufficient to analyze only n pos- 
itive double eigenvalues Q nj k which correspond to k = 
0, . . . , n — 1. The double eigenvalue £l n ,k is semi-simple; 



Now we turn to analysis of the linear stability of the 
modes. Following to the standard procedure, we use the 

substitution q(z,x) = e l/3z [w(x) + u(x)e lu ' z +v*(x)e' lul ' z } the corresponding eigenvectors read p K J k = (w k ,0) T and 
and arrive at the eigenvalue problem 



Lp = up, 



(10) 



where 



L + 2a\w n 



— 2a\w n 



In order to examine splitting of the double eigenvalues, 
we employ Eqs. (|7J|, which yield the following asymptotic 
expansion for the linear stability operator: L = L„ + 
ae 2l Ln^ + o(e 2 ), where 



L = d 2 /dx 2 — b — (x — ia) 2 , and is the Hcrmitian 
adjoint operator. The nonlinear mode w n (x) is unstable 
if there exists an eigenvalue w such that Im uj < 0. 

It is straightforward to check the properties of the op- 
erator L as follows. If u) is an eigenvalue of L with an 
eigenvector (u(x), v (x)) T , then — w* is also an eigenvalue 
with an eigenvector (v*(x), u*(x)) T . Employing the sym- 
metry of the nonlinear modes [w n (x) = iv„(— x)] one 
finds that w* is also an eigenvalue with an eigenvector 
{u*{— x), v*{—x)). Also, u) = is always an eigenvalue of 
the operator L. A corresponding to lu — eigenvector 
reads (w n (x), — (x)) T . 



L (2) 



-b^ ] +2\w n \ 2 



W - 21 



(12) 



Following the standard arguments of the perturbation 
theory for linear operators (22| . in order to explore the 
behavior of a double eigenvalue Q n> k we introduce a 2 x 2 
matrix 



(a) p&.pft-> 



(Pi 



>pL j I*> 



<^ 2) pSU.pS*> \ 



.piT) 



(77) rrrry 



(p (77) -(11), 



I 
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where (a, b) = f (x)a(x)dx for any two column vectors 
a and b. If both the eigenvalues of the matrix M n /. 
are real, then the simple eigenvalues emerging from f2 ni /- 
are real, at least for e > sufficiently small. If such 
a situation takes place for all k — 0, 1, . . . , n — 1, then 
one can state that the nonlinear modes w n (x) belonging 
to the nth family are stable in the linear limit. On the 
other hand, if for some k the matrix M n> fe has a complex 
eigenvalue, then the double eigenvalue V, n ^ gives rise to a 
pair of complex conjugated eigenvalues. This is sufficient 
to conclude that the nonlinear modes of the nth family 
are unstable in the linear limit. For n = no double 
eigenvalues £l n ,k exists. Therefore the lowest family n — 
is always stable in the linear limit. 

Taking into account symmetry of the eigenfunctions 
w n (x), one finds that the entries of the matrix M ni £ have 
the form: 



J w 2 dx 



(M n>fe ) 2)2 = - 
(M„, fe ) 2jl = - 
(M„, fc ) li2 



/ \w n \ 2 w 2 



2n-k 



dx 



/03 



,dx 



u 2n-k L 
J wlwlW2n-kdx 



'2n-k 



dx 



Jwlw k W 2n _ ; 



dx 



Jw 2 k dx 



One also observes that all these entries are real. 

Using the above expressions, the matrices M n ,fc as well 
as their eigenvalues can be found explicitly. One observes 
that for any n and k an expression for the eigenvalues 
of the matrix ~M. n ,k contains a term J P rii k(a), where 
Pn,k( a ) is a polynomial with real coefficients. Such poly- 
nomials are different for different n and k, and are com- 
putable explicitly. Their properties (for n = 1,2... 5) 
are summarized in Table |U 

The terms y/ P n ^{a) represent the only possibility for 
the eigenvalues eventually to have a nonzero imaginary 
part. Respectively, splitting of the double eigenvalue 
il n ,k for a = is determined by the sign of P ra ,fc(0), while 
the behavior of Q nt k in the limit a > 1 is determined 
by the sign of the leading coefficient of the polynomial 
Pn,k(ot). 

Turning now to the Table U the following comments 
can be given: (i) the degree of the polynomial P n ^{oi), 
denoted by obeys the relation D n ^ = 12n — 4fc; 

(ii) more importantly, the leading coefficients of all the 
considered polynomials are positive. It means that for 
any n there exists a critical value a^ 1 * such that for all 
a > a° r the nth family is stable in the linear limit, even 
if this family is unstable in the case of the real harmonic 
potential (i.e. for a = 0). For n = and n = 1 the 
critical values are zero: aff = af — 0. Respectively, 
these families are stable in the linear limit both for the 
case of the real harmonic oscillator (a = 0), and in VT- 
symmetric case for any a. For the next families, n = 2 



n 


k 


D 


s 


s 


positive roots 


1 





12 


+ 


+ 


no positive roots 


2 





24 


+ 


- 


0.05 


2.47 2.54 3.21 3.60 




1 


20 


+ 


+ 




no positive roots 


3 





36 


+ 


- 


0.12 






1 


32 


+ 


- 


0.05 


1.68 1.94 3.18 4.17 




2 


28 


+ 


+ 




no positive roots 


4 





48 


+ 


+ 


0.08 


0.14 3.35 3.40 4.77 4.82 




1 


44 


+ 


- 


0.11 






2 


40 


+ 


- 


0.05 


3.64 4.66 




3 


36 


+ 


+ 




no positive roots 


5 





60 


+ 


+ 


0.12 


0.14 




1 


56 


+ 




0.12 


1.74 2.20 5.14 5.24 




2 


52 


+ 




0.10 


2.41 2.58 




3 


48 


+ 




0.05 


2.84 2.92 4.10 5.10 




4 


44 


+ 


+ 




no positive roots 



and n = 3, the Table U yields 



S cr ) 



3.60, and a 



(cr) 



TABLE I: Properties of the polynomials P n ,k(ot). Here D n< k 
is the degree of a polynomial, S n .k is the sign of the lead- 
ing coefficient, s n> k is the sign of the constant term P„^(0). 
Approximate values of all the positive roots are also reported. 



4.17. This means that being unstable in the linear limit 
for a = 0, the latter families become stable in the linear 
limit for a sufficiently large. The same situation takes 
place to the families n = 4 and n = 5. Moreover we 
conjecture that it also persists for all higher families. 



B. Numerical results 

Passing to the numerical study of the stability (see 
Fig. [T]), we first recall some results known for the real 
harmonic potential, which in our model corresponds to 
a = 0. The nonlinear modes that belong to two lowest 
families (n = and n = 1) are always stable. The fami- 
lies n = 2 and n = 3 are unstable in the linear limit and 
for small and moderate values of U. However both the 
latter families become stable if the nonlinearity is suf- 
ficiently strong (for the stability analysis of the modes 
in strongly nonlinear defocusing medium see Q). In the 
defocusing medium, the value of U, which have to be 
exceeded for the families n = 2 and n = 3 to become 
stable, is large and does not belong to the scope of the 
panel a = of Fig. [TJ 

In the next panel of Fig. [T] we consider the case a 
0. 1 f) For this value of a it follows from Table U that for 
any n = 1, 2, . . . , 5 the nth family is stable in the lin- 
ear limit. Turning to stability of the nonlinear modes of 
arbitrary amplitude, we observe that the lowest family 
n = is stable in the whole the considered region of pa- 
rameters. The same situation takes place for the family 
n = 1 but in defocusing medium only. For a = 1 this 
family loses stability at sufficiently strong nonlinearity. 
The most interesting results, however, are obtained for 
the families n = 2 and n = 3. In contrast to their coun- 
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terparts for the real harmonic oscillator, these families 
are stable in the linear limit. Moreover, these families 
remain stable at least for small and moderate values of 
U. In the defocusing medium, the families n — 2 and 
ri = 3 appeared to be stable in the whole the explored 
region. In the focusing medium, we have found the criti- 
cal values of nonlinearity after which onset of instability 
occurs. It is interesting, that in the certain sense the sit- 
uations for the real oscillator and for the PT-symmetric 
one are opposite: for a = the families n — 2 and n = 3 
are unstable in the linear limit but become stable in fo- 
cusing medium for U sufficiently large. Vice versa, for 
a = 0.15, those families are stable in the linear limit 
but lose stability in focusing medium for large U. At 
this stage we emphasize that only finite range of b and 
U has been considered in our numerics, and in princi- 
ple, the families of nonlinear modes may change stability 
for stronger values of nonlinearity, which have not been 
considered here. 

Next, we have considered 'PT-symmetric harmonic po- 
tentials with stronger imaginary component, a = 1 and 
a = 2. One can deduce from Table Q] that for a = 1 the 
families n = 1, 2, ... 5 are stable in the linear limit, while 
for a = 2 the families n = 1, . . .4 are stable in the linear 
limit and the family n — 5 is unstable. However, both for 
a = 1 and a = 2 all the considered families lose stability 
for relatively small values of U. One observes that the 
larger a, the smaller nonlinearity strength is sufficient for 
the destabilization to occur. 



VI. CONCLUSION 

To conclude we have performed the analysis of the 
structure and the stability of the lowest families of non- 
linear modes in the nonlinear Schrodinger equation with 
the parabolic PT-symmetric potential. We have found 
a number of striking features, not observable for the 
cases of conservative and dissipative parabolic potentials. 
Among these features we emphasize transformation of 
the families bifurcating from the different eigenstates of 
the underlying linear problem to the single family; en- 
hancement of the stability in the linear limit comparing 
to the standard case of the real harmonic oscillator; and 
possibilities of proper choices of the strength of the non- 
conservative part a making unstable for a — nonlinear 
modes to become stable in the PT-symmetric case. 



Acknowledgments 

DAZ was supported by Fundagao para a Ciencia 
e a Tecnologia (FCT) under the grant No. 
SFRH/BPD/64835/2009. VVK was supported by 
the FCT under the grant No. PTDC/FIS/112624/2009. 
The authors acknowledge support by the FCT through 
the grant PEst-OE/FIS/UI0618/2011. 



[1] M. Edwards, K. Burnett, Phys. Rev. A 51, 1382 (1995); 
P. A. Ruprecht, M. J. Holland, K. Burnett, and Mark 
Edwards, Phys. Rev. A 51, 4704 (1995); F. Dalfovo and 
S. Stringari, Phys. Rev. A 53, 2477 (1996); V. I. Yukalov, 
E. P. Yukalova, and V. S. Bagnato, Phys. Rev. A 56, 4845 
(1997). 

[2] M. Kunze, T. Kiipper, V.K. Mezentsev, E.G. Shapiro, 

and S. Turitsyn, Physica D 128, 273 (1999). 
[3] E. G. Shapiro and S. K. Turitsyn, Opt. Lett. 22, 1544 

(1997); S. K. Turitsyn and V. K. Mezentsev, JETP Lett. 

67, 640 (1998); S. K. Turitsyn, T. Schafer, and V. K. 

Mezentsev, Opt. Lett. 23, 1351 (1998). 
[4] Yu. S. Kivshar, T. J. Alexander, S. K. Turitsyn, Phys. 

Lett. A 278, 225 (2001); L. D. Carr, J. N. Kutz, and 

W. P. Reinhardt, Phys. Rev. E 63, 066604 (2001); V. I. 

Yukalov, E. P. Yukalova, and V. S. Bagnato, Phys. Rev. 

A 66, 043602 (2002); R. D' Agosta, B. A. Malomed, and 

C. Presilla, Laser Phys. 12, 37 (2002); V. V. Konotop and 

P. G. Kevrekidis, Phys. Rev. Lett. 91, 230402 (2003); P. 

G. Kevrekidis, V. V. Konotop, A. Rodrigues, and D. J. 

Frantzeskakis, J. Phys. B: At. Mol. Opt. Phys. 38, 1173 

(2005). 

[5] G. L. Alfimov and D. A Zezyulin, Nonlinearity 20, 2075 
(2007). 

[6] D. A. Zezyulin, G. L. Alfimov, V. V. Konotop, and V. 
M. Perez-Gara'a, Phys. Rev. A. 78, 013606 (2008). 

[7] M. P. Coles, D. E. Pelinovsky, and P. G. Kevrekidis, Non- 
linearity 23, 1753 (2010). 



[8] G. Theocharis, P. Schmelcher, P.G. Kevrekidis, D. J. 

Frantzeskakis, Phys. Rev. A. 72, 033614 (2005); D. A. 

Zezyulin, G. L. Alfimov, V. V. Konotop, and V. M. 

Perez-Garcfa, Phys. Rev. A. 76, 013621 (2007); Ruiyu 

Hao, Rongcao Yang, Lu Li, and Guosheng Zhou, Opt. 

Commun. 281, 1256 (2008). 
[9] D. A. Zezyulin, G. L. Alfimov, and V. V. Konotop, Phys. 

Rev. A 81, 013606 (2010). 
[10] see e.g. Special issue of J. Phys. A: Math. Gen. 39 (2006); 

ibidem 41 (2008). 
[11] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 

5243 (1998). 

[12] C. E. Riiter, K. G. Makris, R. El-Ganainy, D. N. 
Christodoulides, M. Segev, and D. Kip D. Nature Phys. 
6, 192 (2010). 

[13] A. Ruschhaupt, F. Delgado, and J. G. Muga, J. Phys. A: 
Math. Gen. 38, L171 (2005); S. Longhi, Phys. Rev. Lett. 
105, 013903 (2010). 

[14] C. M. Bender, D. C. Brody, and H. F. Jones, Phys. Rev. 
D 70, 025001 (2004). 

[15] Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. 
N. Christodoulides, Phys. Rev. Lett. 100, 030402 (2008); 
Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. 
N. Christodoulides, J. Phys. A 41, 244019 (2008); F. Kh. 
Abdullaev, V. V. Konotop, M. Salerno, and A. V. Yulin, 
Phys. Rev. E 82, 056606 (2010); Sean Nixon, Lijuan Ge, 
Jianke Yang, Phys. Rev. A 85, 023822 (2012). 

[16] Sumei Hu, Xuekai Ma, Daquan Lu, Zhenjun Yang, 



7 



Yizhou Zheng, and Wei Hu, Phys. Rev. A 84, 043818 
(2011). 

[17] Zhiwei Shi, Xiujuan Jiang, Xing Zhu, and Huagang Li, 

Phys. Rev. A 84, 053855 (2011). 
[18] Zhenya Yan, Bo Xiong, and W. M. Liu, 

arXiv:1009.4023 ^1 [cond-mat. quant-gas]; V. Achilleos, 

P. G. Kevrekidis, D. J. Frantzeskakis, and R. Carretero- 

Gonzalez, larXiv:1202. 1310V 1 [nlin.PS]. 
[19] X. Zhu, H. Wang, L. X. Zheng, H. Li, and Y. J. He, Opt. 

Lett 36, 2680 (2011). 
[20] Keya Zhou, Zhongyi Guo, Jicheng Wang, and Shutian 

Liu, Opt. Lett. 35, 2928 (2010); H. Wang and J. Wang, 

Opt. Express 19, 4030 (2011); Zhien Lu and Zhi-Ming 



Zhang, Opt. Express 19, 11457 (2011). 
[21] F. Kh. Abdullaev, Y. V. Kartashov, V. V. Konotop, and 

D. A. Zezyulin, Phys. Rev. A 83, 041805(R) (2011); D. 

A. Zezyulin, Y. V. Kartashov, V. V. Konotop, Europhys. 

Lett. 96, 64003 (2011). 
[22] T. Kato, Perturbation Theory for Linear Operators, 

(Springer- Verlag, Berlin, 1966). 
[23] M. Znojil, Phys. Lett. A 259, 220 (1999). 
[24] C. M. Bender and H. F. Jones, J. Phys. A: Math. Theor. 

41 244006 (2008). 
[25] P. M. Morse, H. Feshbach, Methods of theoretical physics. 

Part 7, (New York, McGraw-Hill, 1953). 



